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SECTION  1 


INTRODUCTION 

There  exists  in  the  literature  several  papers  on  the  flight 
mechanics  of  reentry  vehicles.  Unfortunately,  there  is  little 
in  the  literature  that  focuses  on  the  ascent  performance  of  a 
booster  configuration  from  launch  to  apogee.  Consequently, 
there  is  a  need  for  both  approximate  analytical  and  exact 
numerical  trajectory  solutions.  The  former  could  be  used  for 
preliminary  design  purposes.  The  latter  could  be  used  for 
performance  trade-offs  of  finalized  configurations. 

This  investigation  has  two  objectives.  The  first  is  to 
develop  approximate  analytical  trajectory  solutions.  The  second 
is  to  determine  optimal  trajectories  for  a  specified  criterion. 
Numerical  results  will  be  for  a  specific  configuration.  The 
approximate  solutions  are  segments  for  an  ascent  trajectory  that 
can  be  patched  together  to  yield  a  complete  ascent  trajectory. 
The  segments  are  maximum  lift  with  power-on,  constant  flight 
path  with  power-on,  and  zero  angle-of -attack  with  both  power-on 
and  power-off.  For  the  optimal  trajectories,  the  payoffs  are 
maximum  speed  at  burnout  or  apogee  and  minimum  fuel  or  time  to 
reach  apogee. 

In  Section  II,  the  assumptions  employed  in  the  investiga¬ 
tion  are  presented.  Sections  III  and  IV  contain  the  approximate 
and  optimal  trajectory  results. 


SECTION  II 


ASSUMPTIONS 


Spherical  Symmetry 

The  modeling  of  the  atmospheric  density  and  the  accelera¬ 
tion  of  gravity  can  be  greatly  simplified  if  it  can  be  assumed 
that:  1)  the  density  is  a  function  of  only  the  altitude  and  2) 

the  acceleration  of  gravity  is  a  function  of  only  the  radial 
distance  from  the  center  of  the  planet.  The  earth's  surface  is 
an  oblate  spheroid  of  eccentricity  0.00335.  The  ratio  of  the 
radius  at  the  poles  to  that  at  the  Equator  is  0.99665.  The 
difference  in  the  surface  radius  at  the  equator  and  either  pole 
is  approximately  seventy  thousand  feet.  The  ratio  of  the  ac¬ 
celeration  of  gravity  at  the  poles  to  that  at  the  equator  is 
approximately  1.0018.  For  the  purpose  of  this  effort,  it  suf¬ 
fices  to  treat  the  earth  as  spherical 

The  acceleration  varies  inversely  proportional  as  the 
square  of  the  radius  from  the  center  of  earth.  The  nominal 
thickness  of  the  earth's  atmosphere  is  of  the  order  of  275,000 
feet.  The  ratio  of  the  acceleration  of  gravity  at  the  outer 
edge  of  the  atmosphere  to  that  at  the  earth's  surface  is  ap¬ 
proximately  0.972.  In  this  study  the  acceleration  of  gravity 
will  be  assumed  to  be  constant.  The  sea  level  value  of  32.174 
feet/second/second  will  be  used. 

Exponential  Atmosphere 

The  relationship  for  the  rate  of  change  of  density  with 
respect  to  altitude  can  be  approximated  by 

dp 

dh  -  B*  (1) 

where 


The  density  relation  is  exact  across  layers  where  the  tempera¬ 
ture  is  constant.  In  layers  where  the  temperature  is  linear 
with  respect  to  altitude.  Equation  (1)  is  approximate.  If  an 
approximate  value  is  used  across  the  latter  layer,  then  Equation 
(1)  takes  on  the  following  from  for  all  layers. 

B(h  -  hD ) 

_  _  _  -  /  T  \ 


y  =lnp ,  A  =  In  pDe 


k-  "*) 


Then  Equation  (3)  can  be  rewritten  as 

y  =  A  +  Bh  (5) 

Approximation  for  A  and  B  can  be  determined  by  doing  regression 
analysis  on  Equation  (5).  Using  the  1962  U.S.  Standard 
Atmosphere  as  the  basis  gives  the  results  in  Table  1. 

The  accuracy  of  these  approximations  are  presented  in  Table  2 
where  the  error  is  defined  as 


aACT  -  aTHEOR 


%  ERROR  = 


xlOO 


and  a^£,p  is  the  1962  density  ratio  and  is  aPProximate 

value.  The  assumed  form  for  the  density  ratio  for  the  sub¬ 
sequent  performance  analysis  is 

Bh  +  A 


For  high  speed  trajectory  analysis  it  will  be  assumed  that 
the  aerodynamic  coefficients  can  be  predicted  by  Newtonian 
theory.  This  theory  is  based  upon  the  assumption  that  the 
pressure  at  any  point  on  the  windward  side  of  the  vehicle  can  be 
predicted  by 


C  =  K(n*  i )  2  (7) 

P 

Let  the  scaler  product  be  defined  by  sinS  ,then 

C  =  ksin2S  (8) 

P 

The  drag  and  lift  components  due  to  the  pressure  force  are 
evaluated  from 


L 


q  "  J 

C  (k-n) 

1  P 

dA 

(9) 

D 

q  ~  . 

C  (  i  •  n ) 

1  P 

dA 

(10) 

Since 

1  c 

:  cosS 

(11) 

Equations  (9) 


L 

q 

D 

q 


K 

K 


and  (10)  become 

sinzScosSdA 

sin3SdA 


( 12) 

(  13) 


If  the  angle  6  is  constant  everywhere  over  the  vehicle  then  the 
drag  and  lift  coefficients  become 


Ksin  ScosS 


( 14) 


CD  =  qA  =  Ksln  6  1151 

An  additional  component  should  be  added  to  the  drag  coefficient 
to  account  for  bluntness  and  viscous  effects,  thus 


Cp  =  Ksin3S  + 


The  increment  is  a  strong  function  of  nose  bluntness  and 

o 

generally  is  also  altitude  dependent. 

The  approximation  is  going  from  Equations  (12)  and  (13)  to 
Equations  (14)  and  (15)  is  presented  in  the  appendix  where 

S  =  a  +  aQ  (17) 

Seven  different  configurations  were  examined.  The  correlation 
shews  that  Equations  (14)  and  (16)  are  adequate  for  the  purpose 
of  this  investigation. 

Previous  studies  have  shown  that  for  reentry  trajectory 
analysis,  the  lift-to-drag  ratio  is  an  important  parameter.  Fo 
ascent  performance,  the  lift  coefficient  versus  angle-of -attack 
is  an  important  function  for  high  thrust  to  weight  rocket  en¬ 
gines.  Consider  the  aerodynamic  ratio  E  defined  by 


L  Ksin2ScosS 

£  =  —  =  - 

D  Ksin3S  +  CT 


sin2  ScosS 


sin  S  +  C 


where 


The  maximum  value  of  L/D,  E  ,  is  obtained  by  optimizing  with 
respect  to  S.  Differentiating  Equation  (18)  with  respect  to  S 
and  setting  the  resulting  expression  equal  to  zero  gives  the 
necessary  condition  for  E 

sin3SA  +  3Csin2S*  -  2C  =  0  (20: 


This  equation  can  be  easily  solved  for  S  as  a  function  of  C. 

In  general,  however,  we  are  interested  in  determining  the 
parameters  S*  and  C  as  a  function  of  E*.  Solving  Equation  (18) 
for  C  gives 


n  .  2  rA  3C* 

C  =  —*sin  6  cosS  -  sin  6 
E 


(21) 


Substituting  Equation  (21)  for  C  into  Equation  (20)  gives  a 

A  A 

relation  between  E  and  6 


2  -  3sinZ&* 

.  *  7*“ 

3sinS  cosS 


Using  double  angle  formulas,  this  reduces  to 


E*  = 


^  +  cos2S* 
sin26A 


This  is  easily  solved  for  SA  as  a  function  of  EA 


2  _  * 
cos  S 


I  *  2 

'3  +  2E  (E  +  9 
2 ( EA  2  +  1) 


A  A 

Given  E  ,  5  is  determined  from  Equation  (24)  and  C  -vs  obtained 
from  Equation  (21).  These  parameters  and  the  ratio  of  at  5 
to  Cn  are  presented  in  Table  3. 


A 

For  a  parabolic  polar,  CD<S  )/CD  is  2. 

o 

The  maximum  value  of  the  lift  coefficient  is  obtained  from 
optimization  of  Equation  (14)  with  respect  to  S.  The  necessary 
condition  is 

2sinScoszS  -  sin3S  =  0  (25) 

The  solution  for  S  is 

S  =  sin-1 ( 2/3) 1/2  =  54.74°  (26) 

The  value  of  maximum  is 

Cr  =  0.3849K  (27) 

^MAX 

In  Figures  1  through  3,  the  aerodynamic  parameters  are 
presented.  In  Figures  1  and  3,  6  equals  a. 


The  drag  coefficient  can  be  written  in  polar  form  by 
eliminating  S  from  the  CL  and  relations.  Expanding  Equation 
(14)  gives  a  cubic  equation  in  sinzS.  The  solution  of  this 
cubic  equation  is 

1 

sin2S  =  3*(1  -  cosy  +  J3sinp)  (28) 

where 


(  29) 


Substituting  for  sina  in  Equation  (16)  gives 


wV.V.V.' 


Cjj  =  Cp  +  Khj-  (1  -  cosy  +  J3  sinpi) 


For  small  values  of  C^/K,  the  parameters  become  approximately 


sinS  =  JiT 


Cj,  -  C„  +  7 

D  D  „i  /  z 

o  K 


The  optimal  value  for  is 


*  1/2  2/3 

Cr  =  (2K  Cn  ) 

L  Uo 


The  optimal  value  for  CD  is 


CD  =  3CD 


The  maximum  value  of  E  is 


1  f4K 
3  \CZ 


In  Table  4,  exact  aerodynamic  parameters  from  Equations  (14)  and 
(16)  and  approximate  values  from  previous  relations  are 
presented.  It  can  be  seen  that  a  polar  of  the  form  in  Equation 
(32)  gives  fairly  good  agreement  for  E*  when  E*  is  greater  than 
1.5.  If  the  trajectory  requires  very  accurate  data  for  CL  and 
Cp,  then  the  approximate  relations  should  not  be  used.  Accurate 
data  should  be  used  with  accurate  interpolation  techniques 
employed. 


N.'  \ 


Propulsion  Characteristics 


Rocket  propulsion  characteristics  will  be  used  for  the 
performance  analysis.  The  mass  burning  rate  is  determined  from 

WPROP 

W  = t r -  (36) 

CB 

The  burn  time,  tQ,  is  determined  from 

D 

IT0TAL 

tfi  =  -  (37) 

where  I^TAr  is  the  total  impulse  and  has  units  ( force )x(  time) . 
Table  5  presents  propulsion  characteristics  for  several 
boosters.  The  speed  change,  AV,  corresponds  to  the  gravity  free 
solution  in  a  vacuum. 

Equations  of  Motion 


From  Reference  2,  the  equations  of  motion  for  the  trajec¬ 
tory  analysis  are 

dv  1 

=  ~'F,p  -  gsiny  +  w2  rcos0  (  sinycos0  -  cosys  inijjsin<|> )  (38) 

dy  1  V2 

V—  =  -*FjjCOsa  -  gcosy  +  —  -cosy  +  2u>V  costy  cos<J> 

(  39) 

+  w2 rcos<J) ( cosy  cos<$  +  siny  3iniJ>  sin<J>) 
di|>  1  FNslnCT  V2 

Vdt  =  ~  —  •  cosycosi|>tan$  +  2wV(  tanysimjicos^ ) 

(  39) 

w2  r 

-  sin$)  -  — — —  •  costysin4>cos$ 


de 

VcosycosU) 

cTE 

rcos<(> 

(42) 

d<J> 

VcosysiniJ) 

dt 

r 

(43) 

ft 

=  Tcos(a  +  i>v) 

-  D 

(44) 

fn 

=  Tsin(a  +  6v> 

+  L 

(45) 

Some  of  the  terms  in  Equations  (38)  through  (40)  can  be  dropped. 

2 

The  transport  acceleration,  ui  r,  is  less  than  0.4  percent  of  the 
acceleration  of  gravity.  The  Coriolis  acceleration,  2uV,  is 
approximately  eleven  percent  of  the  acceleration  of  gravity  and 
therefore  is  important  for  long  range  high  speed  flight.  For 
low  speed,  short  range  trajectory,  this  term  can  be  neglected. 
For  example,  a  boost  and  coast  to  apogee  trajectory  can  neglect 
the  earth's  rotation.  Also  for  boost  and  coast  to  apogee 
problems,  the  term  Vz/r  in  Equations  (39)  and  (40)  can  be 
neglected  since  these  types  of  trajectories  are  over  3hort 
ranges.  Thus,  the  earth  can  be  treated  as  flat. 

In  the  next  section,  we  turn  our  attention  to  the  boost  and 
coast  to  apogee  problems. 


SECTION  III 


NOMINAL  BOOST  AND  CLIMB  PERFORMANCE 

Both  ground  launches  and  air  launches  are  important 
problems.  The  primary  difference  between  them  is  the  boundary 
conditions.  Ground  launch  initial  conditions  are  zero  speed  and 
ninety  degree  launch  angles.  Air  launches  correspond  to  the 
speed  of  the  launch  platform  and  zero  flight  path  angle.  This 
effort  has  concentrated  on  air  launches  only. 

A  typical  trajectory  consists  of  the  separation  from  the 
carrier,  propulsion  ignition,  pitch  up  to  a  specified  attitude, 
and  then  an  acceleration  to  burnout.  This  may  be  followed  by 
the  separation  of  the  booster  and  then  a  second  boost  to  a 
second  burnout.  After  the  boost  phase  or  phases,  the  vehicle 
coasts  and  decelerates  to  apogee.  An  impulse  may  follow  to  give 
the  desired  orbital  speed. 

Five  different  boost  and  climb  segments  are  of  interest. 


are 

as  follows: 

a . 

maximum  lift,  power-on 

b. 

constant  flight  path. 

power-on 

c . 

zero  angle-of -attack. 

power-on 

d. 

zero  angle-of -attack. 

power -of  f 

e . 

optimum  boost  climb 

The  last  problem  is  addressed  in  Section  IV.  This  problem 
is  difficult  and  merits  attention  in  a  section  by  itself.  The 
difficulty  is  determining  the  optimal  controls.  The  first  four 


make  up  a  typical  sequence  from  launch  to  apogee.  Each  of  the 
first  four  problems  is  solved  in  the  following  subsections. 


Flight  is  in  the  equational  plane  and  the  earth's  rotation 
is  neglected.  The  equation  for  the  rate  of  change  of  r,  r,  is 
replaced  with  the  time  rate  of  change  of  the  density.  The 
result  is 

ct  =  BaVsin^  (46) 

where  B  is  defined  in  Equation  (6).  The  remaining  differential 
equations  are 

T 

W  =  -  7 —  (47) 

iSP 

1 

V  =  -*Ft  -  gsinv  (48) 

.  1  9 

*  =  ^’FN  '  V'COSY  (49) 

V 

0  =  —’cosy  (50) 

where 

Frj,  =  Tcos(a  +  Sv)  -  D  (51) 

Fn  =  Tsin( a  +  6  )  +  L  (52) 

Nondimensional  transformation  are  introduced  in  order  to 
normalize  the  variables  to  the  order  of  one.  Let 

u  =  V/V( o ) ,  u  =  W/W( o ) ,  0  =  t/t_  (53) 

K 

We  are  free  to  select  the  reference  time,  tR"  The  differential 
equations  become 


(54) 


a  =  Xx<ju  siny 


w 


u 


T 

A  (--cosS  -  d  -  siny) 

3  U) 


Y  =  - (~*sinS  +  n  -  cosy) 

U  U) 

0  =  X^UCOSY 

where 


A  =  tRBV(0) 
fcRT 

X2  =  W(0)lsp 


A 
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V(0) 

V(0)tR 


R 


WPROP 
W(  0) 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


S  =  a  +  S 

v 

D  L  T 

d  =  W'  n  =  M'  T  W(o) 


(63) 

(64  ) 


In  Equation  (50),  r  has  been  approximated  by  R.  In  Equation 
(62),  the  parameter  A4  is  very  small.  Thus  the  change  in  the 
downrange  is  small.  In  Equations  (54)  through  (58),  differen¬ 
tiation  is  with  respect  to  the  nondimensional  time  0. 


The  approach  is  to  develop  approximate  closed  form  solu¬ 
tions  for  the  different  trajectory  segments.  These  approximate 
solutions  are  compared  with  the  exact  solutions  obtained  by 
numerical  integration  of  the  differential  equations  of  motion. 


This  comparison  offers  some  evidence  as  to  the  fidelity  of  the 
approximate  solutions.  As  expected,  some  accuracy  is  lost 
through  introducing  assumptions  that  make  it  possible  to 
analytically  integrate  the  differential  equations.  But,  one 
advantage  is  that  an  approximate  solution  can  be  easily  obtained 
through  substitution  of  specific  values  for  the  parameters.  We 
turn  our  attention  to  the  first  problem.  The  booster  charac¬ 
teristics  correspond  to  configuration  17  in  Table  5. 

Maximum  Lift,  Power-On 

This  segment  corresponds  to  the  initial  rotation  after 
launch.  The  initial  conditions  are  the  launch  conditions.  The 
terminal  conditions  are  free.  Integrating  Equation  (55)  gives: 


w  =  1  -  (65) 

Assume  that  the  flight  path  angle  is  small  and  the  angle-of- 
attack  is  constant.  Also,  assume  that  the  change  in  the  drag 
coefficient  is  negligible,  which  is  a  questionable  assumption, 
and  the  change  in  the  density  is  negligible,  which  is 
reasonable.  Then  the  drag  ratio  in  Equation  (64)  is  ap¬ 
proximately 

u2 

d  =  — -d(0)  (66) 

w 


Equation  (56)  is  approximately 

x 

u  =  — •  tcosS  -  u2d(0) 


(67) 


Combining  with  Equation  (55)  and  integrating  gives 


where 


s  = 


TCOS  6 


"duT  '  x6  =  2  ~ •J'rd<0)c°afi 

2 


In  a  manner  like  that  for  d,  n  is  approximately 


u 

n  =  n( o) - 

u> 


(6 


(7 


At  high  angles-of -attack,  the  right  side  of  Equation  (57)  is 
approximately  n,  since 

T 

-•sinS  -cosy  (7 

approximate  forms  for  Equations  (56)  and  (57)  gives 

n(  0 )  u 

- -  (7 

tcos 6  -  d(0)u 

Integrating  gives 

>  =  y<0>  -  2  E(0)lrV  (0)  l7- 


n  >  > 

Combining 

dy 

du 


where 


E(C) 


n( 0 )  L( 0 )  CL(0) 

dToT  =  dToT  =  CD(0) '  fT  =  TCOs6  “  d(0)u 


(7- 


Combining  Equations  (54)  and  (56)  gives  approximately 


x  ouy 
1 


tcosS  -  d ( 0 ) i 


If  the  change  in  the  weight  is  neglected,  i.e.,  assume  w  =  1, 
then  Equation  (75)  can  be  integrated 


a(  0 ) 


1  Xi  T  i  1  fT 

U0rHV(0)  +7E(0)lnft(0)Jlnj-J51 


E(  0 )  l 


8d( 0 )  x 


( Hr] 2  -  [ 


-  In  fT(i 


Combining  Equations  (58)  and  (56)  gives  approximately 


du  x. 


xcosS  -  d( 


0)u2j 


Assuming  w  =  1  permits  integration 


6  0(O)  2d( 0 ) *x3'ln  fT(0) 


In  Table  6,  the  comparison  between  the  exact  and  ap¬ 
proximate  solutions  is  presented  for  the  first  ten  seconds.  The 
flight  path  reached  approximately  thirty  degrees.  In  light  of 
the  fact  that  the  exact  solution  included  acceleration  through 
the  transonic  region,  the  agreement  is  quite  good.  The  density 
used  in  the  results  correspond  to  tabular  values  rather  than 
Equation  (6).  Equation  (6)  gives  a  lower  density,  approximately 
5.5  percent.  Using  this  approximate  value  in  the  equation  for 
d(o)  results  in  little  difference  for  the  speed,  altitude,  and 
downrange.  But  the  flight  path  angle  is  lower  by  approximately 


5.5  percent.  The  drag  and  lift  coefficients  are  the  initial 
values . 

Constant  Flight  Path,  Power-On 

This  is  a  path  that  could  be  followed  once  the  initial 
rotation  is  completed.  The  objective  for  this  type  of  trajec¬ 
tory  is  to  establish  an  attitude  or  flight  path  angle  that  will 
result  in  attaining  specified  apogee  conditions.  Generally,  the 
angle-of -atcack  is  small  since  high  longitudinal  accelerations 
are  desired  for  constant  flight  path  angles.  Equation  (57) 
gives 


-rsinS  +  un  =  wcosy(O) 


(79) 


The  right  side  decreases  with  increasing  time;  therefore,  the 
angle-of -attack  must  decrease  as  the  speed  increases.  The 
initial  conditions  are  those  at  the  beginning  of  the  segment 
which  are  the  final  conditions  at  the  end  of  the  previous  seg¬ 
ment,  e.g.,  the  initial  rotation. 

It  is  assumed  that  the  change  in  the  drag  in  Equation  (56) 
is  negligible.  Combining  Equations  (55)  and  (56)  and  integrat¬ 
ing  gives 

X  3  X 

u  =  1  -  —  •  FxcosS  -  d(0)lln  w  +  —  •<(*>  -  l)siny(O)  (80) 

z  *-  -*  2 


Combining  Equations  (54)  and  (55),  substituting  Equation  (80), 
and  then  integrating  gives 


a 


In 


o(0) 


1  /  X3 

-•  sinv(  0 ) | (w  -  1)  -  —  •  j^tcosS  -  d(G)j 


(81) 


A. 

r  -|  13 

|wlnu)  -  (uj  -  1)J  +  —  •  —  •(u  -  DsinyfO) 


Combining  Equation  s(58)  and  (54)  and  integrating  gives 


0 


0(0) 


+ 


x 

— • coty( 0 ) • In 

A 


o 

o(0) 


(82) 


The  exact  and  approximate  solutions  are  presented  in  Table 
7.  The  agreement  is  surprisingly  good  in  light  of  the  assump¬ 
tion  made  about  the  drag. 

Zero  Angle-of -Attack.  Power-On 

Along  this  path,  the  lift  is  zero  and  the  drag  coefficient 

is  minimum  since  CD  equals  .  This  corresponds  to  maximum 

o 

acceleration,  at  least  on  a  point  performance  basis.  Over  the 
total  segment,  however,  the  integration  of  the  drag  may  be 
larger  relative  to  positive  angles-of -attack.  This  is  at¬ 
tributed  to  the  fact  that  a  zero  angle-of -attack  trajectory  is 
lower  in  the  atmosphere;  thus,  the  dynamic  pressure  may  be 
higher . 


It  will  be  assumed  that  the  change  in  CD  is  negligible  and 

o 

the  flight  path  angle  is  small.  The  flight  path  angle  decreases 
as  the  vehicle  accelerates.  It  is  assumed  that  the  change  in 
the  drag  is  negligible.  The  speed  increases,  but  the  density 
ratio  decreases.  Equation  (56)  reduces  to 


UVl»ft  H.i  1.4  l 


Combining  with  Equation  (55)  and  integrating  gives 


u  =  1 


x  -  d ( 0 ) I lnu 


The  inverse  of  this  is  required  in  order  to  integrate  y 


u  =  e 


C(  1  -  u) 


where 


a  t  -  d  ( 


Combining  the  approximate  differential  forms  for  u  and  y,  sub 
stituting  for  u,  and  then  integrating  gives 


oo  (-Cu)1  -  (-C)1 


Y  =  y  (  0)  -  — -Ce  lnu  +  E 
z  i  =  l 


i  -  i ! 


Combining  Equations  i 54 )  and  (56),  substituting  Equation  (85) 
and  integrating  gives 


a  t'  C  r 

3  1  CT  e  -Cu  -C 

— - In— 777  - e  LU(Cu  +  l)y  -  e  (C  +  l)y(O) 

A  x  CT(0)  c2  L 


2C  r 
e  i 


-H 


e'2Cu-e  2e)  +  lnu  +  I 


oo  (-1)1  (  2C  ) 1 


i  •  i  ! 


(u  -  1) 


where 


Combining  the  approximate  differential  equations  for  6  and  u  and 
integrating  gives 


x  e 

4 


0  =  0(0)  + 


A3fTC 


(C  +  1 ) e-C  -  (Cu  +  1 ) e 


-Cul 


(90) 


The  solution  for  the  nondimensional  weight  is 


UJ  =  1  -  Xz0  (91) 

The  exact  and  approximate  solutions  are  presented  in  Table 
8.  The  agreement  is  quite  good. 


Zero  Angle-of -Attack.  Power-Off 


This  segment  is  the  path  after  burnout  up  to  apogee.  The 
objective  is  to  reach  apogee  with  minimum  loss  of  energy.  This 
corresponds  to  minimum  drag  which  is  zero  angle-of -attack  on  a 
point  performance  basis.  Assume  that  the  aerodynamic  forces  are 
negligible  relative  to  the  gravity.  Equations  (56)  and  (57)  can 
be  combined  to  give 

du 

t~  =  utanv  (92) 


This  is  easily  integrated  to  give 

cos^( 0 ) 

u  =  - 

cosy 


(93) 


Substituting  for  u  in  Equation  (57)  and  integrating  gives 


tanv  -  tan-y(O) 


s*3 

cos-/(  0 ) 


(94) 


Combining  Equations  (54)  and  (56)  and  integrating  gives 


In 


a 

o(0) 


>95) 


Combining  Equations  (58)  and  (57)  and  integrating  gives 


A 

0  =  0(0)  -  — •  cos2y(  0 )  j^tany  -  tany(0)J  (96) 


The  comparison  between  the  exact  and  approximate  solutions 
is  presented  in  Table  9.  The  agreement  is  surprisingly  good 
considering  the  drag  is  neglected.  It  is  desirable,  however, 
that  the  drag  be  negligible  for  the  power-off  coast  to  apogee. 

The  exact  solution  reaches  apogee  at  135.6  sec. 

Complete  Trajectories 

With  the  relations  derived  earlier,  it  is  easy  to  develop  a 
complete  trajectory  from  launch  to  apogee.  The  initial  condi¬ 
tions  are  the  launch  conditions.  Thereafter,  the  initial 
conditions  are  the  end  conditions  from  the  previous  segment. 

For  comparison  purposes,  an  approximate  trajectory  of  130 
seconds  was  developed  and  compared  with  the  exact.  The  ap¬ 
proximate  solutions  were  based  upon  the  methods  developed  in  the 
previous  sections.  The  first  ten  seconds  of  the  trajectory  was 
maximum  lift.  The  next  twenty  seconds  was  constant  flight  path 
angle  with  power-on.  The  next  thirty  seconds  was  zero  angle-of- 
attack  with  power-on.  The  last  seventy  seconds  was  zero  angle- 
of -attack  with  power-off.  The  exact  and  approximate  solutions 
are  presented  in  Table  10. 

Next  we  turn  our  attention  to  optimal  boost  climb  trajec¬ 
tory  analysis. 


SECTION  IV 


OPTIMAL  BOOST  CLIMB  PERFORMANCE 

The  objective  of  this  analysis  is  to  determine  the 
aerodynamic  controls  as  a  function  of  time  such  that  a  boost 
ascension  trajectory  occurs  in  some  optimal  manner.  In  all 
problems,  the  angle-of -attack  is  a  control  variable.  In  some 
problems,  the  thrust  vector  angle  is  also  a  control  variable. 
When  it  is  not,  the  thrust  vector  angle  will  be  set  equal  to 
zero . 


The  optimal  controls  are  a  function  of  the  performance 
payoff  criterion.  There  are  three  payoffs  that  are  of  interest: 
the  first  is  maximum  speed,  the  second  is  minimum  propellant  to 
reach  specified  terminal  conditions,  and  the  third  is  minimum 
time  to  reach  apogee.  The  second  payoff  is  equivalent  to  maxi¬ 
mizing  the  burn  out  weight  since  the  weight  varies  linearly  with 
time . 


The  relevant  differential  equations  are  Equations  (54) 
through  (57).  Equation  (58)  can  be  dropped  since  6  is  not  a 
payoff  criterion;  it  does  not  contain  the  control,  and  the  other 
differential  equations  do  not  contain  9.  Since  x2  is  a  con¬ 
stant,  Equation  (55)  can  be  replaced  with  its  integral 


UB0 

0  i  P  1  —  =  P 
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»B0  ^  *> 


'  97) 


where  tDr.  is  the  burnout  time.  It  is  assumed  that  the  thrust 
and  aerodynamic  coefficients  are  independent  of  altitude  and 
Mach  number.  This  appears  to  be  a  gross  approximation,  but 
generally  the  acceleration  from  subsonic  to  hypersonic  speeds  is 


very  rapid  for  high  thrust  to  weight  conf igurations .  At  hyper¬ 
sonic  speeds,  the  aerodynamic  coefficients  will  be  a  function  of 
only  the  angle-of -attack. 


The  boost  is  at  maximum  thrust  until  burnout.  For  some 
problems,  this  is  the  terminal  condition.  For  other  problems, 
the  vehicle  coasts  to  apogee  after  burnout.  This  corresponds  to 

\Of)  =  0  (98) 

The  optimal  trajectory  problems  are  identified  in  Table  11. 
The  solutions  to  Problems  4  and  9  are  trivial.  If  the  final 
conditions  are  free,  then  minimum  propellant  or  minimum  time  is 
zero.  Thus,  the  initial  and  final  times  are  the  same.  Since 
the  propellant  used  is  a  function  of  time.  Problems  5,  6,  7,  and 
8  are  equivalent  to  10,  11,  12,  and  13,  respectively. 

The  following  examples  will  be  addressed  later  on  in  this 
section.  In  Problem  1,  the  payoff  is  the  final  speed  and  there 
are  no  constraints.  The  remaining  problems  are  apogee  problems 
and  are  defined  as  follows  in  Table  12. 


Aerodynamic  Relations 


Since  the  drag  and  lift  coefficients  are  respectively  even 
and  odd  functions,  rewrite  Equations  (14)  and  16)  as  follows 

CD  =  K | sin3a |  +  CD  (99) 

o 

=  K  |  sina.  |  sinacosa  (100) 

where  |x|  means  absolute  value  of  x.  For  the  configuration  to 

be  used  in  this  section,  C„  =  1.035  and  K  =  35.30.  In  order  to 

Do 

develop  necessary  and  sufficient  conditions  for  the  optimal 
controls,  we  need  first  and  second  derivatives  of  and  with 
respect  to  a..  The  derivatives  are 


2  2 

C'l  =  K(2cos  a  -  sin  a)|sina| 


2  2 

=  K(2cos  a  -  7sin  <x)cosa  a  >  0 


=  -K(2cos2a  -  7sinza ) cosa  a  <  0 


C  =  3C 
a 


CD  3CL 
oca  a 


(101) 


(102) 


(103) 


( 104) 


The  coefficients  and  their  derivatives  are  presented  in  Figures 

4  through  7.  The  range  of  a  is  that  from  -Cr  to  +  Cr 

^MAX  LMAX 

The  optimal  control  problem  is  considered  next. 

Optimal  Control  Problem 

The  approach  to  determining  the  solution  of  the  optimal 
control  problem  is  given  in  Reference  2.  The  optimal  control 
problem  is  the  maximization  of  the  Hamiltonian  with  respect  to 
the  control  variables.  The  Hamiltonian  is  defined  as  follows 


H  =  Pu  +  P  y  +  P  a 
u  Y  a 


( 105) 


where  P  ,  P  ,  P  are  the  costate  variables  associated  with  the 
u  y  o 

state  variables  u,  y,  o.  The  necessary  condition  for  a  costate 
variable  x  (u,  y,  or  o)  on  an  optimal  trajectory  is 


(  106) 


Consequently,  three  additional  differential  equations  result. 

If  the  thrust  vector  angle  is  free,  there  are  three  control 
variables.  Since  6  is  the  sum  of  a  ar.d  only  two  can  be 


Vi*  *  / 


•*  J 

v  v  *  ■ 


considered  as  independent.  The  choice  is  <x  and  S  .  If  the 

v 

thrust  vector  is  fixed,  Sv  =  0,  then  only  a  is  free.  Necessary 
conditions  for  both  cases  will  be  addressed.  Consider  first 
free  Sv  and  a.  Introduce  two  control  functions 


F(6)  =  CrptcosS  +  P^sinS) 

VI  K 


G( a)  =  F(a)  -  CD  +  PR  CL 


( 107 ) 


(  108) 


where 


S  =  a  +  §  ,  C_  =  —  ,  P_  -  — =- 
v  T  qs  R  uP 


The  part  of  the  Hamiltonian  that  is  a  function  of  Sv  is 


H<6  )  =  X  P  •  — -F(S  ) 

V  3  U  W  V 


The  part  of  the  Hamiltonian  that  is  a  function  of  oc  is 


( 109) 


( 111) 


H(a)  =  X  P  •^n’G(ot) 
3  u  W 


If  S  is  unconstrained,  optimal  5  ,  S  ,  is  the  solution  of 
v  r  v  v 

3F(  6  ) 
v 


(  111 ) 


(112) 


Substitution  gives 


tani  a  +  S  )  =  Pr, 

V  K 


(113) 


If  a  is  unconstrained,  optimal  <x,  <x  ,  is  the  solution  of 


3G(ou 


(  114  ) 


Since  Fft  =  =0  according  to  Equation  (112),  then  a  must  be 

the  solution  of 


CD  +  PR  CL  “  0 
a  a 


(116) 


Substituting  Equations  (103),  (100),  and  (101)  gives 


-3tana*  +  PR(2  -  tanza*)  =  0 


(117) 


We  can  solve  for  oc  as  a  function  of  PR  or  PR  as  a  function  of 
a*.  We  choose  the  latter  so  that  we  can  examine  the  Hamiltonian 
as  a  function  of  a*  and  S* 


3  tana 

R  2  -  tanza* 


(118) 


Thus,  once  a  is  determined,  P0  is  determined  from  Equation 

A  ** 

(118)  and  S„  is  obtained  from  Equation  (113).  In  Figure  8,  PR 

V  ^  A  ^ 

is  presented  as  a  function  of  a  .  In  Figure  9,  S  is  presented 
as  a  function  of  a  .  The  functions  F  and  G,  evaluated  at  a  and 
6*,  are  presented  in  Figures  10  and  11  as  functions  of  a  and 


In  order  that  the  Hamiltonian  be  a  local  maximum,  it  is 


necessary  that 


Both  relations  are  evaluated  for  the  optimal  controls.  In 
Figures  12  and  13,  Fg  g  and  are  presented.  Since  both 

partial  derivatives  are  negative,  it  follows  from  Equations 
(110)  and  (111)  that  along  the  optimal  trajectory 

P  >  0  (121) 

u  — 

From  Equation  (109)  and  Figure  8,  it  follows  that  if  <x  <.  0, 

PD  <  0,  therefore,  P  <0.  If  a  >  0,  PD  >  0,  therefore,  P  >  0. 
Furthermore,  if  a  =  0,  P^  =  0. 

If  the  thrust  vector  angle  is  fixed,  assume  =  0,  then 
along  the  optimal  trajectory 


G  =  F  -  Cn  +  PD  Cr  =0  (122) 

a  a  D  R  L 

a,  a 

In  this  case,  Fft  *  0.  The  optimal  angle-of -attack  is  the  solu¬ 
tion  of 


CT(  -sinoi*  +  PRcosaA)  -  +  Pp  c£  =  0 

a  a 


( 123) 


Substituting  Equations  (103),  (100),  and  (101)  gives 


A  A  .  A  A  A 

C,p  (-sina  +  PpCosa  )  -  3K|sina,  |sma  cosa 


+  PRK(2cos2aA  -  sin2<xA  )  |  sincx*  |  =  0 


( 124) 


The  optimal  angle-of -attack  is  a  function  of  PR,  K,  and  C,p. 
Previously  it  was  a  function  of  only  Pp.  By  substituting  tan¬ 
gent  trigonometric  relations  for  sine  and  cosine  terms.  Equation 

•  A 

(124)  can  be  transformed  to  a  cubic  equation  in  tana 


Cp( -x  +  PR)(1  +  x2>  "  3K | x  j  x  +  PRK(2  "  x2 )  ! x |  -  0  (125) 


where 


Equation  (125)  can  be  solved  for  a  as  a  function  PQ,  K,  and  C-, 
or  PR  solved  in  terms  of  a  ,  K,  and  C^,.  Thus 


C,jOc(l  +  x2)  +  3K|x|x 


CT(1  +  x2)  +  K( 2  -  x2) | x | 


( 127) 


PR  is  presented  in  Figure  14.  In  Figure  15,  G(a  )  is  presented. 

H  A 

Gq ux(ot  )  is  presented  in  Figure  16. 

For  maximum  Hamiltonian,  Relation  (119)  must  hold.  From 
Equation  (111)  it  follows  that  along  the  optimal  trajectory 


P  >  0 
u  - 


(128) 


This  agrees  with  the  previous  case.  Thus  for  both  situations, 
free  or  fixed  thrust  vector  angle.  Equations  (121)  and  (128) 
hold. 

We  next  turn  our  attention  to  the  development  of  the  op¬ 
timal  trajectory  solutions. 

Optimal  Trajectory  Solutions 

There  are  two  ways  to  determine  the  optimal  trajectory. 

The  first  way  involves  determining  the  costate  differential 
equations  according  to  Equation  (106)  and  then  integrating  these 
equations  between  appropriate  boundary  conditions  The  second 
way  is  an  extension  of  the  first  and  involves  developing  dif¬ 
ferential  equations  for  the  controls.  These  equations  are  also 
a  function  of  Equation  (106). 


The  first  approach  will  be  employed  here.  Two  differential 
equations  that  involve  the  costate  variables  are  required.  Only 
two  are  required  since  the  costate  variables  can  be  normalized 


with  respect  to  P  .  The  differential  equations  are  in  terms  of 

u 

PR  and  P  where 


P 

a 


(129) 


From  Equation  (106),  the  differential  equations  are 


-•A3cosy  -  PgX^COSY  -  --X3Pr^.cosS  +  dj 

1  2  fT  ) 

-  -'PR  x3  I  w*  sinS-  n  -  cosy  I  +  PRPsXiasinY 


(130) 


S  1 


x  d  —  -  -  -  — • x  nPn 

3  la  u  )  u  3  R 


psl  ‘  XiP S usinY 


(131) 


"u'X3  PSPR'W-sinS  +  u*X3PSPRCOSy  +  PSXiasin> 

These  equations  and  the  trajectory  equations  are  integrated 
forward  in  time  to  the  burnout  condition.  Optimal  controls  are 
imposed  everywhere  along  the  trajectory.  Initial  guesses  for 
Pg(0)  and  PR(0)  or  a(0)  are  required.  For  problem  1,  the  final 
flight  path  angle  and  altitude  are  free,  thus 


W  =  0 


pR(tf)  =  0 


( 132) 


Results  are  presented  in  Figures  17  through  23.  The  differences 
for  free  Sv  or  Sv  =  0  are  negligible  for  the  speed  and  altitude. 
The  payoffs  for  Problem  1  are  presented  in  Table  13.  Since  the 
difference  in  the  speed  for  the  first  problem  is  negligible,  the 
remaining  problems  will  be  constrained  by  Sv  =  0. 

For  problems  2,  3,  and  4,  the  optimal  trajectory  solutions 
are  presented  in  Figures  24  through  27.  The  initial  values  for 
a*  and  Pg  are  presented  in  Table  14  along  with  the  constraints 
and  optimal  payoff.  For  reference  purposes,  the  maximum  speed 
for  Problem  1  and  6V  =  0  is  V(t^)  =  6583  ft/sec  and  =  60 


seconds.  The  maximum  speeds  for  Problems  2  and  3  are  less  than 
Problem  1.  Also,  the  speeds  for  Problem  3  are  less  than  Problem 
2.  The  minimum  time  solutions  for  Problems  4  are  less  than  the 
corresponding  maximum  speeds  for  Problem  3.  The  flight  path 
angle  solutions  show  that  the  minimum  time  solution  pitches  up 
to  a  higher  flight  path  angle.  Once  this  is  reached,  the  op¬ 
timal  angle-of -attack  solution  shows  that  there  is  a  rapid  pitch 

over  until  the  angle-of -attack  for  -Cr  is  reached.  This 

^MAX 

angle-of -attack  is  followed  until  apogee  is  reached.  This  is 

one  of  the  primary  differences  between  the  controls  in  the  two 

problems.  For  minimum  time  to  apogee,  the  initial  pitch  up  is 

raster  and  the  pitch  down  is  also  rapid  until  -Cr  is  reached. 

^MAX 

Recall  that  the  angle-of -attack  was  limited  to  angles  between 

-Cr  and  +Cr  .  A  second  difference  is  that  the  minimum 
MAX  ^MAX 

time  trajectory  may  reach  apogee  before  burnout,  whereas,  the 
maximum  speed  solution  at  apogee  never  occurs  before  burnout;  it 
must  occur  at  or  later  than  burnout. 


SECTION  V 


CONCLUSIONS 

Nominal  and  optimal  boost  climb  trajectories  were  developed 
for  an  air  launched  vehicle.  The  approximate  nominal  solutions 
were  in  good  agreement  with  the  solutions  obtained  by  numerical 
integration  of  the  differential  equations  of  motion.  The  ad¬ 
vantage  of  the  approximate  solutions  is  a  rapid  estimate  of  the 
ascent  of  a  rocket  powered  vehicle. 

For  the  optimal  tra jectories ,  three  payoffs  were  con¬ 
sidered;  maximum  speed  at  burnout  or  apogee,  minimum  fuel  to 
apogee,  and  minimum  time  to  apogee.  The  last  two  problems  are 
equivalent  since  the  fuel  burning  rate  i3  constant.  Thrust 
vectoring  was  considered,  but  the  difference  relative  to  no 
vectoring  was  negligible  for  the  problem  of  maximizing  the  speed 
it  burnout.  As  a  result,  all  subsequent  problems  assumed  no 
thrust  vectoring. 

In  general,  the  optimal  angle-of -attack  decreased  with 
increasing  time.  The  change  was  greater  for  the  minimum  time 
solution.  For  this  case,  a  fast  pitch  up  occurred  initially  and 
then  the  angle- of -attack  rapidly  decreased  to  the  angle-of- 

attack  for  -Cr  .  This  gave  the  fastest  pitch  over  to  apogee. 

MAX 

The  minimum  time  solution  may  reach  apogee  before  burnout. 
The  emphasis  is  on  reaching  apogee  before  burning  all  propel¬ 
lant.  For  the  problem  of  maximum  speed  at  apogee,  all 
propellant  is  used  prior  to  or  at  apogee.  The  reason  for  this 
is  that  a3  long  as  there  is  propellant,  the  vehicle  will  con¬ 
tinue  to  accelerate. 

The  analysis  has  been  limited  to  constant  thrust  and  burn¬ 
ing  rate  and  a  single  booster.  If  the  constants  are  variable, 
then  relationships  as  a  function  of  the  appropriate  variables 


must  be  included.  This  will  change  the  costate  differential 

equations  if  the  engine  variables  are  dependent  upon  any  of  the 
state  variables.  If  there  is  more  than  one  booster,  there  may 
be  coast  period  between  boosts.  Care  must  be  exercised  to 
insure  that  this  possibility  is  properly  accounted  for. 
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APPENDIX 


CORRELATION  BETWEEN  EXPERIMENTAL  AND 
THEORETICAL  AERODYNAMIC  DATA 

Seven  different  configurations  were  examined.  Maximum 
lift-to-drag  ratio  varied  from  approximately  one  to  over  four. 
The  seven  configurations  are  the  X24C,  Reference  3,  Super 
Orbital  Re-entry  Test  Integrated  Environment  (SORTIE),  Reference 
4,  Shuttle  Orbiter  Vehicle  102,  Reference  5,  Advanced  Manned 
Interceptor  (AMI),  Reference  6,  X20  Dynasoar,  Reference  7, 
Aerothermodynamic/Elastic  Structural  Systems  Environmental  Tests 
(ASSET),  Reference  8,  and  an  optimized  configuration  from 
Reference  9.  The  configurations  are  illustrated  in  Figure  A-l 
where  the  AMI  has  the  highest  experimental  maximum  L/D  ratio  and 
the  SORTIE  is  the  lowest.  In  all  but  the  optimized  configura¬ 
tion,  either  CL  and  CQ,  or  CN  and  CA  were  available.  If  the 
latter  were  available,  then  the  data  were  transformed  to  CL  and 
CQ.  For  the  optimized  configurations,  only  L/D  distributions 
were  available. 

The  experimental  data  were  obtained  from  wind  tunnel  tests. 
The  theoretical  distribution  were  assumed  to  correspond  to  the 
relations  in  Section  II,  Equations  (14)  and  (16) 


CL  =  Ksinz5cos6 

(  14) 

CD  =  Ksin3S  +  CD 

o 

( 16) 

will  be  assumed  that 

S  =  a  +  a 

o 

(  17) 

There  are  different  ways  for  correlating  the  data.  The 

approach  here  was  to  determine  K,  CD_  ,  and  <xQ  such  that  a  sat is - 

o 

factory  fit  to  L/D  and  C^  was  obtained.  For  each  data  set,  a 


subset  of  the  data  was  selected.  The  error  is  defined  as  fol¬ 
lows 


ERROR 


( theor ) -E( exp) 


( theor ) -C^ 


( exp) 


( A-l ) 


sinzScosS 
sin3S  +  C 


E( exp) 


KsinzScos6-CL( exp) 


where  is  the  first  point  and  N2  is  the  last  point.  The 
objective  was  to  minimize  the  error  with  respect  to  the 
parameters  C  and  K.  The  parameter  C  was  obtained  from  the 
minimization  of  the  squared  error  between  the  therotical  and 
experimental  values  for  L/D  and  wa3  accomplished  by  employing  a 
golden  section  search  technique.  The  minimization  of  the  error 
with  respect  to  K  was  obtained  by  minimizing  the  error  in  CL> 
For  the  optimized  configuration  only  C  was  obtained.  For  all 
but  this  configuration,  <xQ  was  assumed  to  be  zero.  For  this 
configuration  it  appeared  that  aQ  =  4°  gave  the  best  fit  to  the 
daca. 


The  results  are  presented  in  Figures  A-2  through  A-8.  The 
emphasis  was  on  a  satisfactory  fit  to  the  L/D  data  in  the  neigh¬ 
borhood  of  E*.  Consequently,  some  of  the  data  does  not  fit  well 
at  larger  angles-of -attack.  If  maneuvers  at  large  angles-of- 
attack  are  required,  than  the  emphasis  should  be  on  the  region 
of  interest.  It  can  be  seen  that  the  experimental  and  theoreti¬ 
cal  distributions  match  quite  well  in  the  neighborhood  of  E*. 

The  values  of  N  (total  number  of  experimental  points),  N. , 

A  A  ■** 

N7,  Cn  ,  K,  6  ,  and  E  are  presented  in  Table  A.  The  solutions 


for  6  and  E  were  obtained  from  Equations  (20)  and  (23).  For 
the  optimum  vehicle,  it  was  assumed  that  E*  =  3,  olq  =  4°,  and  6 * 
and  C  were  computed  from  Equations  (24)  and  (21).  for  the 
purpose  of  this  investigation,  the  theoretical  distributions  for 
CL  and  Cq  are  satisfactory. 


Table  1 


ATMOSPHERIC  PARAMETERS 


A  B 


0  <  h  1 

25,000  feet 

0.9858  x  10~2 

-0.3209 

X 

10 

25,000  < 

:  h  <  150,000  feet 

0.4231 

-0.4695 

X 

10 

150,000 

±  h  <  250,000  feet 

-0.6150 

-0.3956 

X 

10 

250,000 

1  h 

3.734 

-0.5704 

X 

10 

Table  2 

EXPONENTIAL  DENSITY  ERRORS 


Altitude 

Error 

(ft) 

(%> 

0 

-1.0 

25,000 

-1.0 

50,000 

4.1 

75,000 

-3.7 

100,000 

-1.3 

125,000 

-3.4 

150,000 

-3.0 

175,000 

-0.8 

200,000 

3.9 

225,000 

5.8 

250,000 

-4.4 

275,000 

1.8 

300,000 

-3.1 

Table  3 


AERODYNAMIC  PARAMETERS 


ea 

S  * ( Deg ) 

CD  /K 

cd(6*)/ce 

1 

29.32 

.0916 

2.28 

2 

17.57 

.0159 

2.73 

3 

12.24 

.0051 

2.86 

4 

9.34 

.0021 

2.92 

EXACT  AND  . 

Table  4 

APPROXIMATE  AERODYNAMIC  PARAMETERS 

CD  /K 
o 

EXACT 

E*  C*/K 

CD/CD 

0 

APPROXIMATE 

E*  C*/K 

cA 

LD 

0.3675 

0.5 

0.320 

1.74 

0.74 

0.814 

0.0916 

1.0 

0.209 

2.28 

1.17 

0.323 

0.0342 

1.5 

0. 132 

2.57 

1.63 

0.167 

0.0159 

2.0 

0.087 

2.73 

2.10 

0.100 

0.0051 

3.0 

0.044 

2.86 

3.07 

0.047 

0.0022 

4.0 

0.026 

2.92 

4.05 

0.027 

*  P 


Table  5 


SAMPLE  BOOSTER  PROPULSION  CHARACTERISTICS 


Booster 

T/W( 0 ) 

tg( sec ) 

AV(f t/sec) 

WPR0P/W( 0) 

1 

2.58 

38.8 

4200 

0.44 

2 

2.94 

3.9 

5300 

0.54 

3 

15.69 

5.4 

3500 

0.40 

4 

18.40 

3.2 

2400 

0.30 

5 

1.87 

28.9 

2100 

0.31 

6 

3.75 

26.9 

4300 

0.45 

7 

20.17 

4.3 

3600 

0.40 

8 

9.42 

5.5 

1900 

0.25 

9 

21.41 

3.4 

2800 

0.31 

10 

16.40 

3.0 

1800 

0.20 

11 

2.18 

60.0 

5900 

0.49 

12 

2.84 

60.0 

8900 

0.64 

13 

2.66 

59.3 

9900 

0.68 

14 

3.00 

60.1 

9500 

0.67 

15 

2.54 

60.1 

7200 

0.56 

16 

2.71 

59.0 

7700 

0.58 

17 

2.51 

60.0 

8200 

0.69 

Table  7 


CONSTANT  FLIGHT  PATH,  POWER-ON  COMPARISONS 

The  conditions  from  the  exact  solution  at  the  end  of  the  initial 
pitch-up,  t  =  10  seconds,  are  the  beginning  of  the  constant  flight 
path,  power-on  trajectory. 


O  0 


V(0) 

=  1248 

ft/sec,  yc 

=  29. 

25  , 

h(  0  )  =  42,500 

ft,  0(0) 

=  .026 

W(  0) 

=  13,272  lb,  W  = 

172.5 

lb/3ec,  T  =  37,600 

lb,  tR  = 

50  sec 

S  =  0, 

cD(°) 

=  0. 

95,  S  =  10.9 

ft2 

EXACT 

APPROXIMATE 

At 

V 

h 

0 

V 

h 

0 

0 

1248 

42,500 

.026 

1248 

42,500 

.026 

10 

1903 

50,300 

.063 

1958 

50,300 

.064 

20 

2690 

61,800 

.117 

2809 

61,900 

.121 

30 

3701 

77,600 

.192 

3855 

78,000 

.200 

40 

5089 

99,400 

.295 

5190 

100,000 

.307 

50 

7095 

129,500 

.437 

7044 

130,000 

.453 

Table  8 


ZERO  ANGLE-OF-ATTACK,  POWER-ON  COMPARISONS 

The  conditions  from  the  exact  solution  at  t  =  30  seconds  are  the 
initial  conditions  for  the  zero  angle-of -attack,  power-on  trajec¬ 
tory. 


V(0)  =  2707  ft/sec,  y(0)  =  30°,  h(0)  =  62,000  ft,  0(0)  =  0.118°, 


W(0) 

=  9787 

.5  lb. 

W  =  172.5 

lb/sec , 

CD(0) 

=  0.75, 

S  =  10.9 

ft2. 

T 

=  38,500 

lb,  a  - 

°,  tR 

=  30  3ec 

EXACT 

APPROXIMATE 

At 

V 

X 

h 

0 

V 

X 

h 

0 

0 

2707 

30.0 

62,000 

0.118 

2707 

30.0 

62,000 

0.118 

5 

3193 

27.3 

69,100 

0.154 

3268 

26.9 

69,400 

0.159 

10 

3753 

24.9 

76,700 

0.196 

3886 

24.3 

77,300 

0.208 

15 

4405 

23.0 

84,900 

0.247 

4573 

22.1 

85,900 

0.266 

20 

5173 

21.2 

93,900 

0.307 

5349 

20.3 

95,000 

0.333 

25 

6087 

19.8 

103,700 

0.379 

6238 

18.7 

104,800 

0.413 

30 

7193 

18.6 

114,600 

0.464 

7281 

17.3 

115,400 

0.504 

Table  9 


ZERO  ANGLE-OF-ATTACK  POWER-OFF  COMPARISONS 


The  initial  conditions  correspond  to  the  exact  solution  at  bur¬ 
nout,  t  =  60  seconds. 


0  O 


V(0) 

=  7193 

■  ft/sec 

,  y(0)  = 

18.6  , 

h(  0 )  = 

117,400 

ft,  0(0)  = 

0.464 

W  = 

4612  pounds,  a 

=  0,  tR 

=  70  sec. 

EXACT 

APPROXIMATE 

At 

V 

X 

h 

0 

V 

X 

h 

0 

0 

7193 

18.6 

114,600 

0.464 

7193 

18.6 

114,600 

0.464 

10 

6990 

16.3 

135,800 

0.648 

7097 

16.1 

135,900 

0.651 

20 

6863 

13.9 

153,800 

0.830 

7014 

13.6 

154,100 

0.837 

30 

6770 

11.5 

168,900 

1.010 

6946 

11.0 

169,000 

1.023 

40 

6699 

9.0 

180,900 

1.190 

6891 

8.4 

180,600 

1.210 

50 

6646 

6.5 

189,900 

1.369 

6852 

5.7 

189,100 

1.397 

60 

6609 

4.0 

196,000 

1.548 

6827 

3.1 

194,300 

1.583 

70 

6588 

1.4 

199,100 

1.727 

6817 

0.4 

196,400 

1.770 

The  exact  solution  reaches  apogee  at  135.6  sec. 

V  =  6582  ft/sec 
h  =  199,600  ft 
0  =  1.827° 

The  approximate  solution  reaches  apogee  at  t  =  131.3  sec. 

V  =  6817  ft/sec 
h  =  196,400  ft 
0  =  1.795° 


Table  10 


COMPLETE  TRAJECTORY  COMPARISONS 


EXACT 

APPROXIMATE 

t 

v 

X 

h 

0 

V 

X 

h 

0 

0 

700 

0 

40 , 000 

0 

700 

0 

40,000 

0 

10 

1248 

29.2 

42,500 

0.026 

1311 

29.5 

42,800 

0.030 

30 

2690 

29.2 

61,800 

0.117 

2862 

29.5 

62,900 

0.127 

60 

7193 

18.6 

114,600 

0.464 

7376 

17.3 

117,500 

0.524 

130 

6580 

1.4 

199,100 

1.727 

7043 

-0.5 

192,200 

1.873 

OPTIMAL  CONTROL  PROBLEMS 


MAXIMUM 

PAYOFF 

CONSTRAINTS 

V(tf ) 

none 

V(tf ) 

It 

4J 

>■ 

0 

V(tf ) 

ytt^.)  =  0,  h(tf) 

specif ied 

W(tf ) 

none 

W(  tf ) 

y(tf)  = 

0 

W(tf ) 

y(tf)  =  0,  h(tf> 

specified 

W<  tf ) 

y(tf)  =  0,  V(tf) 

specified 

W<  tf ) 

y(tf>  =  0,  V(tf)  & 

h(tf)  specified 

-fcf 

none 

y(tf)  = 

0 

-fcf 

y(tf)  =  0,  h(tf) 

specif ied 

-fcf 

y(tf)  =  0,  V(tf) 

specif ied 

y(tf)  =  0,  V(tf)  & 

h(tf-)  specified 

Table  12 


SELECTED  OPTIMAL  CONTROL  PROBLEMS 
y(  t^. )  =  0 


PROBLEM 

PAYOFF 

CONSTRAINT 

2 

V(tf ) 

none 

3 

V(tf ) 

h(tf ) 

4 

-tf 

h(  tf ) 

Table  13 

PROBLEM 

1  RESULTS 

Sv 

A 

V  (tf) 

a*<0) 

Pg  ( o  > 

free 


6609  ft/sec 


31.9° 


-5.80 


Table  14 


OPTIMAL  TRAJECTORY  RESULTS 


Problem 

<x*(0) 

ps(0) 

Payoff 

tf ( SEC ) 

Constraint 

2 

35.86 

-6.601 

V(tf >=6256ft/sec 

60 

none 

3a 

36.05 

-6.555 

=6235f t/sec 

60.4 

h(tf)=100, 000ft 

3b 

35.95 

-6.093 

=5935f t/ sec 

99.9 

=150, 000ft 

4a 

40.20 

-7.732 

t £=48. 7sec 

— 

=100, 000ft 

4b 

39.30 

-6.331 

=74 . 3sec 

— 

=150, 000ft 

Table  A 

AERODYNAMIC  PARAMETERS 

Vehicle 

N 

N1  N2 

C°o 

K 

SA(Deg) 

ea 

AMI 

5 

3  3 

.0081 

4.279 

8.86 

4.23 

OPTIMUM 

6 

- 

- 

- 

12.25 

3.00 

X24C 

7 

4  5 

.0222 

2.705 

14.29 

2.54 

DYNASOAR 

7 

3  3 

.0488 

2.411 

19.06 

1.83 

SHUTTLE 

11 

1  11 

.0543 

2.821 

18.77 

1.86 

ASSET 

6 

3  3 

.  1276 

3.303 

23.43 

1.43 

SORTIE 

7 

5  7 

.1181 

1.581 

29.11 

1.09 

OPTIMUM  VEHICLE 


MI 


X24C 


*1 


DYNASOAR 


"7 


>- 


■  4 


SHUTTLE 


ASSET 


^  \ 

,  '  SORTIE 
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